Early Universe Tomography with CMB and Gravitational Waves 
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We discuss how one can reconstruct the thermal history of the Universe by combining cosmic 
microwave background (CMB) measurements and gravitational wave (GW) direct detection exper- 
iments. Assuming various expansion eras to take place after the inflationary reheating and before 
Big-Bang Nucleosynthesis (BBN), we show how measurements of the GW spectrum can be used to 
break the degeneracies associated with CMB data, the latter being sensitive to the total amount of 
cosmic expansion only. In this context, we argue that the expected constraints from future CMB 
and GW experiments can probe a scenario in which there exists late-time entropy production in 
addition to the standard reheating. We show that, for some cases, combining data from future CMB 
and GW direct detection experiments allows the determination of the reheating temperature, the 
amount of entropy produced and the temperature at which the standard radiation era started. 



PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



Our understanding of the evolution of the Universe is 
now becoming clearer owing to precise cosmological ob- 
servations such as cosmic microwave background, large 
scale structure (LSS), type la supernovae and others. 
From such observations, we can obtain information about 
the current energy budget and the history of the Uni- 
verse. In particular, the evolution after the time of Big- 
Bang Nucleosynthesis to the present is relatively well un- 
derstood. On the other hand, one can also probe the evo- 
lution during inflation since cosmic density fluctuations, 
which can be probed by CMB and LSS, are considered 
to be initially generated during that epoch. 

Compared to the the above mentioned eras, the evolu- 
tion, or thermal history of the Universe during the period 
after inflation to BBN is relatively unexplored, certainly 
due to the lack of associated cosmological observables. 
Although, in the standard scenario, the Universe is con- 
sidered to be radiation-dominated until BBN (precisely 
speaking, until the radiation-matter equality epoch) after 
the inflaton reheating, the thermal history can be more 
complicated. In theories beyond the standard model of 
particle physics such as in supersymmetric models and 
string theory, there can exist some scalar fields (other 
than the inflaton, as for instance moduli field) that are 
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long-lived and can dominate the energy density of the 
Universe. Their decay may also produce huge amount of 
entropy thereby influencing the early universe history. 

In the light of these considerations, it would be worth 
investigating how one could probe the thermal history 
during these epochs. In fact, some authors have investi- 
gated this issue by using observations of CMB [H-Q and 
direct detection of gravitational waves (GW) in Refs. 0- 
0]. From CMB observations, we can probe the amplitude 
of the primordial scalar and tensor fluctuations as well 
as their scale dependencies around the so-called pivot 
scale. Notice that the time at which such a reference 
scale exited the Hubble radius during inflation depends 
on the amount of cosmic expansion from Hubble exit 
to the present times, which of course includes all of the 
above-mentioned post-inflationary eras. As a result, by 
measuring the primordial power spectra in a given infla- 
tionary model, we can obtain information on the amount 
of the total cosmic expansion, i.e. the integrated ther- 
mal history since the end of inflation. This can be also 
applicable in GW direct detection experiments through 
measurements of the amplitude and scale dependence of 
tensor fluctuations. Not only that, direct detection of 
GWs could be used to probe the background evolution 
as the shape of the GW's spectrum is very sensitive to 
it. Thus, in the inflationary framework, detection/non- 
detection of GWs can give invaluable information on the 
thermal history of the Universe thereby allowing a "to- 
mography" of these eras. 

In this paper, we investigate this issue by complement- 
ing observations of the CMB and GWs, paying particular 
attention to the period from the end of inflation to BBN. 
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For this purpose, we first recall how one can constrain the 
thermal history of the Universe from these experiments. 
Although we have not detected any gravitational waves 
yet, CMB observations, such as those from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) |8|, are pre- 
cise enough to already give some constraints within some 
inflationary models [3|]. However, in the near future, a 
direct detection of GWs could be achieved for some in- 
flationary models that would allow to combine both CMB 
and GWs experiments. To see this in an explicit man- 
ner, we investigate the projected constraints on the ther- 
mal history of the Universe from future observations of 
CMB such as CMBpol and the PLANCK satellite com- 
bined with future direct detection GW experiments such 
as BBO @, DECIGO [Io| and Ultimate-DECIGO 

The organization of this paper is as follows. In the 
next section, we give a brief description on how CMB 
and GWs can probe or constrain the thermal history after 
inflation and justify their complementarity. We also give 
the current constraints on the thermal history within the 
so-called large field model of inflation coming from CMB 
using WMAP data. Then in Section ITO1 we present our 
forecasts derived from a Fisher matrix analysis based on 
the above-mentioned future CMB and GWs experiments. 
We conclude in the last section. 



states h\ in Fourier space as 

dk 3 

(2^f 



h ij (T),x)= j J0j3j2 h ^ k '- 



X ikx 



(3) 



where e^- are the polarization tensors satisfying e^e l y — 

At first order in the slow-roll formalism, t he p rimordial 
power spectrum for the scalars is given by |12j | 
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while the tensor spectrum (sum of polarization included) 
reads 



2Hj 
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In these equations, M 2 j = 1/(8ttG) stands for the re- 
duced Planck mass, C = — 2 + In 2 + -fE — 0.73 with ^e 
being the Eulcr constant and e; are the slow-roll param- 
eters which are defined as 1 



II. CMB AND GWS AS A PROBE OF THE 
THERMAL HISTORY OF THE UNIVERSE 



In this section, we describe how one can probe the 
thermal history of the Universe with CMB and GW ob- 
servations in the context of inflationary cosmology. 

We assume that inflation is the origin of both 
scalar and tensor perturbations around the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric 



ds 2 = -a 2 (l + 2$)d77 2 



[(l - 2*)<y i; 



hij] dx*dx J 



(1) 

where \& and $ are the Bardeen potential and hij is the 
transverse and traceless spin two fluctuations. If inflation 
is driven by a slowly-rolling scalar field <fi, the quantum 
fluctuations of the field-metric system generate an al- 
most scale invariant power spectrum for both kinds of 
perturbation. The observable quantities are those which 
are conserved on super-Hubble scales, that is the comov- 
ing curvature for scalar perturbations, which reads in the 
longitudinal gauge, 



6<j>(ri,x) 



(2) 



where % = aH is the conformal Hubble parameter and 
a prime denotes derivatives with respect to the confor- 
mal time r]. The tensor modes hij are themselves gauge 
invariant and conserved on super-Hubble scales. It is 
convenient to decompose them on their two polarization 
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where N — In a is the number of e- folds. The spectral 
index defined as n s — 1 = d In Vq / d In k is given by 



l-2e x , 



£2* 



(9) 



The tensor-to-scalar ratio r, which is usually used to 
quantify the amplitude of the tensor mode, is given by 



Vh i« 
vl = 16ei *' 



(10) 



An asterisk "*" indicates that the quantities have to be 
evaluated at the time when the pivot mode fc* crossed 
the Hubble radius during inflation, i.e. the solution of 



a(r)*)H(r)*). 



(11) 



1 In the literature, the slow-roll parameters defined using the po- 
tential for the inflaton V(rp) are a lso used, which are given by 
e v = (l/2)M2j (V'/V) 2 and jjy = M^(V" /V) with a prime 
representing the derivative with respect to (/>. The correspon- 
dence between {ey , r)y} and {ei, £2} at leading order in slow-roll 
parameters are : 



ei = ey, 62 = 2ey — 2r]y. 



(6) 



With these parameters, the spectral index and the tensor- to- 
scalar ratio are respectively given by 



ns = 1 — 6ey„ + 2r)y„, r = 16ey*. 



(7) 
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Here, we neglect the running of the spectral index, which 
will not affect our results for CMB. However, we note that 
such truncation can give large deviation for the power 
spectrum from the one obtained exact numerical calcu- 
lation in some models 1131 . 



A. CMB constraints on the post-inflationary 
universe history 

1. Standard scenario 

The power spectrum functional forms of Eqs. (j4]) and 
([5]) are usually compared to the current CMB data 
to constrain the slow-roll (Hubble flow) parameters, or 
equivalently the spectral index and tensor-to-scalar ra- 
tio [3, H, [T3 - [T9| . This is done by choosing a pivot scale 
k„ in the observable range, typically fc* = 0.05 Mpc -1 . 
However, if one assumes an inflationary model, there is 
much more to say. Indeed, the tensor-to-scalar ratio, 
the spectral index and all other observable quantities are 
completely determined by the inflaton potential V(<p). As 
discussed earlier, they have to be evaluated at the time 
ry*, which can be determined by solving Eq. ([lip. In or- 
der to obtain 77*, it is compulsory to make assumptions 
on the subsequent thermal history of the universe, i.e. 
including at least the reheating and preheating stages. 
In terms of the number of e-folds during inflation, the 
physical pivot wavenumber is given by 



a 



a 



(1 + z end )e N ° 



-N 



(12) 



where z cnd is the redshift at which inflation ended, after 
-/Vend e-folds. As shown in Refs. [lH3|, |20J, a convenient 
way to calculate z en d is to introduce the so-called "re- 
heating parameter" 



-Rrad = 



&cnd / Pcnd 



&reh \ Prch 



1/4 



(13) 



The quantity i? rad encodes all of our ignorance of the sub- 
sequent thermal evolution after the end of inflation and 
quantifies any deviations from a pure radiation era. In 
fact, by assuming instantaneous transition between the 
inflationary epoch to inflaton oscillating era and inflaton 
oscillating to radiation-dominated eras, one has 



1 + ^nd = 



1 



1? 



rad 



Pcnd 
P70 



1/4 



(14) 



where p cn d is the energy density of the universe at the 
end of inflation and p 7o is the energy density of radiation 
today, eventually rescaled by any change in the number 
of gravitating relativistic degrees of freedom. In terms of 
the cosmological parameters today, 



P 7o — QrchPfo = 3Q r 



,h M p 2 ro 



(15) 



where we have defined 



_ <?reh / ffso 
telrch = 



90 \5sroh 



4/3 



(16) 



Here, g B and g respectively denotes the number of en- 
tropic and energetic relativistic degrees of freedom at the 
epoch of interest. Ho and il ro are the Hubble parameter 
and radiation density parameter today. 

As shown in Ref. [2fj| . using energy conservation, 
Eq. (fl~3"|) can be recast into two other strictly equivalent 
forms 



lni? rad = ~(A re h - A cnd )(3u> rc h - 1) 

1 — 3W rcn / /j re h 
— In 



12(l + wT reh ) 



Pcnd 



(17) 



with Wreh standing for the mean equation of state pa- 
rameter during the inflaton oscillating era. Using i? rad , 
Eq. (dl|) is solved for the e-fold time A A* = A* - A ond , 
verifying [20| 



1 



Hi 



AA* = -lni? rad + A --ln. 

4 \M^e u 



4 Vei* K 3-ei cnd 



where the constant An stands for 



An = In 



fc*/qp 
v P70 , 



(18) 



(19) 



Let us emphasize that the right hand side of Eq. (p~5|) 
usually depends on AA* itself, but in a completely alge- 
braic way once the model, i.e. V(<fi), is specified. It can 



also be simplified further by using ei„ 
and Eq. (0} as, 



< 1, ei 



end 



1 



AN. 



lni? rad +A --ln(87r 2 'P, 



1 r 1 V end 

As an example, we have plotted in Fig. Q] the pre- 
dicted values for ei», £2*, as well as the spectral in- 
dex n s and tensor-to-scalar ratio r for the chaotic in- 
flation models with V((f>) oc (f> p . We show the cases of 
p = 1,2,3,4,6 and 10. For each p, there is a range of 
possible AJV, since the reheating should occur from the 
end of inflation to the BBN epoch, whose values are in- 
dicated in the figure at both end points and for each 
case of p. Along with the theoretical predictions, we also 
show la and 2a confidence intervals associated with the 
WMAP7 datajl [H [H[ and Hubble Space Telescope 
(HST) data [23] . From the figure, one can easily see 
that inflation models with p > 3 are excluded by current 
data. Even for the cases with p = 2 and 1, there is a 
lower bound on |AA»| to be consistent with WMAP7, 
which can be translated into the constraints on the ther- 
mal history of the Universe. 



4 




FIG. 1: WMAP7 constraints in the plane (ei„, £2*) (top) and 
(n s ,r) (bottom) compared with the large field model predic- 
tions obtained by solving Eq. (|18|) for various monomial po- 
tentials V(4>) oc (j> p . The annotated values are those of |AiV,,| 
and they range for a reheating occurring as low as BBN to an 
instantaneous reheating after inflation. 



2. Non-standard scenarios 

Up to here, we have assumed the "standard" scenario 
in which once the reheating from the inflaton is com- 
pleted, the Universe becomes radiation dominated until 
the radiation- matter equality z cq ~ 10 4 . However, this 
standard scenario could be modified. For instance, if one 
assumes that, inserted after the reheating era there is a 
phase of evolution dominated by a gravitating source X, 
characterized by an equation of state parameter w x . As 
shown in Ref . |2G | , one can define a parameter i? x exactly 
as in Eq. §T5§ by 



&xend \ Pxcnd 



1/4 



(21) 



"end — > xini", "reh — > xend". From this definition, one 
can check that all equations are unchanged, and in par- 
ticular Eq. ([TBI) , by replacing i? ra d with i? ra di? x . As- 
suming another F-era to take place just after the X-era 
and before the radiation-dominated era, we would reach 
exactly the same conclusions by replacing R la dRx with 
Ri&dRxRy ■ In other words, CMB can only constrain the 
overall thermal history and only feels those parameters, 
Rr&dRxRy ■ • • , multiplied. Let us also notice that the 
correction coefficient entering Eq. (j!5[) is now given by 
Q y end instead of Q ron . 

As a well motivated example, and the one we will be 
discussing in Sec. lIIIl such a situation is typical of scenar- 
ios in which a late-decaying massive scalar field, denoted 
as a hereafter, produce a large amount of entropy well 
after the inflaton reheating. In that case, the X-era is a 
short radiation-dominated era standing just after inflaton 
reheating and before the a domination. For such a sce- 
nario, one has i? x = 1 whereas, the Y-era would precisely 
correspond to the field domination era having w y = 
such that R y can only take negative values (quadratic 
potential). In Fig. [2] to Fig. Uwe have represented the 
WMAP7 constraints on the combination i? ra di? x i? y (for 
i? x = 1) for various large field models, either massive 
as the scenario we are interested in, or for any values 
of p, the power law exponent of the inflaton potential. 
The method we have used is the same as in Ref. Q and 
we do not repeat the details here. As already shown 
in Fig. [TJ large values of |AiV*|, which corresponds to 
RmdRy < 0, give bigger r and more red-tilted spectral 
index n s . Hence smaller values of i? ra d-R y are disfavored, 
which means that the current observations already give 
some constraints on the thermal history of the Universe. 
The posteriors for p en d are also depicted and since it is 
essentially determined by the amplitude of the curvature 
perturbation, it is well bounded [3J. 

In some literature, when one considers a late-time en- 
tropy production scenario, the parameter F is adopted 
to quantify the amount of entropy production, instead of 
R y , and is defined as 



F : 



Syend Qyend 



SyiniQyini 



(22) 



The subscript "yini" and "yend" indicate that the quanti- 
ties are the ones evaluated at the time when K-era starts 
and ends. Contrary to the definition of the i? ra d and R x 
parameters which only require that total energy density 
is conserved, the definition of F is convenient if thermal- 
ization is achieved. In our scenario, at the beginning of 
the Y-era, and also just after its end, the Universe is as- 
sumed to be radiation dominated, and if thermalized, the 
entropy density is dominated by relativistic species. In 
that situation, it is straightforward to show that 

1/4 



for which it is immediate to verify that Eq. (TT71) also ap- 
plies using the mean value w x and by the replacement 



Ry = F- 1 ^ 

where Q is defined as in Eq. 
est. 



Q 



yim 



(23) 



Qycnd 

161) for the epochs of inter 
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FIG. 2: Marginalized posterior probability distributions (solid line) and mean likelihood (dotted line) from WMAP7 and HST 
data for massive inflation with p — 2. The right figure shows the one- and two-sigma confidence intervals of the two-dimensional 
marginalized posterior in the plane (R Ta ,aR y , fc 4 p C nd)- 






FIG. 3: Marginalized posterior probability distributions (solid line) and mean likelihood (dotted line) from WMAP7 and HST 
data for the large field potential with a power law exponent p being varied in 0.2 < p < 5. 



B. Stochastic gravitational waves background 

As discussed in the previous section, CMB can con- 
strain the amount of entropy production but this will 
remain completely degenerated with reheating from the 
inflaton as the only quantities appearing in the deter- 
mination of AN* is the product of R parameters such 
as Rr&dRxRy As we show below, this is not the same 
for stochastic gravitational waves of inflationary origin: 
they feel these parameters in a different way which can be 
used to break the degeneracies, thereby performing the 
tomography of the history of the Universe. In particu- 
lar, direct detection experiments as BBO and DECIGO, 
which probe the frequency range / ~ 0(1) Hz, would 
give new and complementary information with respect 
to CMB. 

In order to discuss the amplitude of stochastic GW, one 



usually uses the spectrum of the energy density of GW 
normalized by the critical energy density p CT n- From the 
pseudo stress-tensor, assuming a stochastic background 
in which spatial and time averages are identical, one 
gets 13, [25 1 



1 dp 



gw 



Pent d In k 



1 

12 
1 

12 



k 
k 

oH 



k 3 . 



A 



(24) 



V 



(obs) 



(k). 



The observed power spectrum "p£ obs ) at the time of in- 
terest (e.g. today) is given by Eq. (|5|) times a transfer 
function encoding the evolution of the sub-Hubble modes 



(obs) 



(k)=V h (k)T 2 (k). 



(25) 
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FIG. 4: Two-dimensional posteriors with la and 2a confidence regions from WMAP7 and HST data in the (p, R Ta dRy) plane 
(left) and in the (7? ra d-R y , K 4 p C nd) plane (right). 



The transfer function T can be evaluated analytically 
since, after inflation, the tensor perturbations are decou- 
pled from other sources and their equation of evolution 
in Fourier space is given by 



h" + 2Uti + k 2 h = . 



(26) 



where for simplicity the polarization index has been 
dropped. Assuming for the moment that the background 
is dominated by a fluid having a constant equation of 
state parameter w, the solution reads 

h(r),k) oc (/by^^/^+^Js^/by). (27) 

For super-Hubble modes (krj -C 1) we recover that h(rj, k) 
is constant while for sub-Hubble wavenumbers one has 



\h(ri, k)\ 2 oc a 

fe»7>l 



(28) 



Since h stays constant on super-Hubble scales, the trans- 
fer function is determined only by the era at which a 
given mode k reenters the Hubble radius. Every mode 
k will then be damped till the present time by a factor 
(afc/ao) 2 where ak is solution of 

k = a k H{a k ) , (29) 

and the scale factor today will be taken as unity clq = 1. 

1. Radiation and matter eras 

We start with the case for which the mode reen- 
ters during radiation era. Assuming an instantaneous 



transition, the result can be approximated by squaring 
Eq. (|29l) and using the first Friedmann-Lemaitre equa- 



tion, H 2 = p/(3M 2 {). One gets 



3M p 2 ,fc 2 



3M p 2 ,fc 2 4 



cq 



where p oc a 4 is used. The solution reads 



ak 



-cq 

eq k 



1 



with 



1 



eq V3M pl 1 



/Pcq 



(30) 



(31) 



(32) 



Z p j J. -T i C q yj fi rQ 

In the rightmost side, f2 mo and ft To are the matter and 
radiation density parameters today and p eq ~ 2p l0 /a^ q is 
the energy density at equality. For all modes entering the 
Hubble radius in the radiation era, Eq. (I31|) immediately 
gives the transfer function for k > fc eq [26j: 



T la d{k) 



rr2 ^eq 



k 



(33) 



Similarly, starting from Eq. ()29[) . the matter era solu- 
tion reads 



T mat (k) =s ^S^e(fc < k cq ) , 



(34) 



where Q(x) is the Heaviside step function. In fact, con- 
sidering a mixture of radiation and matter, one would 
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find 



Teq(fc) 



2k 2 




(35) 



which gives back Eqs. (|3"3"f and in the appropriate 
limits. Although not specified, the above transfer func- 
tions are unity for supcr-Hubble modes, i.e. k < H^ 1 . 
Combined with Eq. (|24[) we recover the well-known re- 
sult [25| that f2 gw is constant for k > k cq and decays as 
k~ 2 (x f~ 2 for modes entering during matter domination. 

The same line of reasoning can be applied to the post- 
inflationary universe assuming the non-standard history. 
Just after reheating we assume the universe to be in a X- 
era (radiation-dominated in our scenario) , then becomes 
F-dominated, e.g. driven by an oscillating scalar field, 
which finally decays into radiation. 



2. Non-standard post-inflationary eras 

If an observable mode today entered the Hubble radius 
during the F-era, Eq. (|29l) can be dealt exactly as the 
radiation-era case, provided w y remains constant. The 
equivalent of Eq. (1501) now reads 



4 



^k 2 y 

3+3u! y k 
Pyend« yen d 



whose solution can be recast into 



^yend 



^yend 



2/(l+3iu y ) 



(36) 



(37) 



This expression can be further simplified by remarking 
that 

1/4 



aoQ^a , (38) 

\PyendJ 



since the end of the F-era matches with the beginning of 
the standard radiation dominated era. The wavenumber 
fcyend correspond to a mode entering the Hubble radius 
just at the end of the F-era, i.e. at the beginning of the 
radiation era: 



Eyelid 



^yend 



pi 



VPy 



1 



\/3M 



-(P7oPycnd) 1/4 , (39) 



pi 



where the last equality comes from Eq. (|38|) . At this 
point, Eq. (|37|) shows that the quantities k ycn d and w y 
are observable and completely determined by the mea- 
surement of f2 gw (fc). As a result, we should try to ex- 
press all quantities in terms of them, and in particular 
the redshift at which the F-era ended. From Eq. (|38|). 
one gets 



1 



*yend 



= V3 



fcycnd-Mpl 

-1/2 
Pfo 



(40) 



V3fc ycnd M p i V k ) 
+ Q(k < k ycad ) . 



We finally get the transfer function during the F-era 

T y {k) ~ ^ - (*>=*) e(k> fcycnd ) 

(41) 

It remains now to deal with the modes entering the 
Hubble radius before, i.e. either during X domination or 
during reheating. The calculations are again the same 
although the scale factor at the end of the X-era reads 



_ ^xend t-> 

^xend — Oyend — ^ycnd-^Cy 

Oycnd 



/ \ !/ 4 

/ Pycnd \ 

V Pyini / 



(42) 



where Eq. (|2ip has been used for a xe nd/%end- One can 
further simplify this expression by using Eq. (|38j) to get 



^xend — -^y 



Pxcnd 



1/4 



(43) 



Defining /c xcn d the wavenumber of a mode entering the 
Hubble radius at the end of the X-era, we have 



, ^xend , 

™xcnd — TtFTZ V Pxcnd 



i?v 



pl 



V3M, 



pi 



■(p 70 p xc „d) 1/4 , (44) 



such that the corresponding redshift can be expressed in 
terms of observable quantities as 



1 + z xcnd = V3 



R y Pla 



(45) 



The transfer function during this era is again given by 
\ak\ and reads 



Tx(fc) 



R 2 -1/2 

-f«,yP70 



Pl 



V3 fcxendM] 
+ 0(fc < /Cxcnd) 



2/(l+3ui x ) 

1 6(/c > fc xcnd ) 



(46) 

Finally, in order to deal with modes entering the Hub- 
ble radius during reheating (inflaton oscillation domi- 
nated era), we similarly express a rc h in terms of a xcn d 
to get 



^rch R^Ry 



Pjyo 
Preh 



1/4 



(47) 



Again, the wavenumber crossing the end of reheating is 
given by 



, _ a rch , _ 



RxRy ,~ \l/4 



V3M pl v ^ V3M pl 
and the redshift at which reheating ended reads 



1 + z Ich = V3 



RlR y p 70 



(48) 



(49) 
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Therefore, the transfer function during reheating is 

+ 9(fc < fc rch ) , 

which shows the influence of i? x and i? y , making them 
observable in the gravitational wave spectrum. However, 
contrary to CMB, they are no longer degenerated with 
i? ra d. Conversely, using both CMB and 0, gw , one expects 
to be able to disambiguate the effects of i? ra d, R* and 
R y . The same reasoning could be extended to another 
Z-era inserted somewhere. 

To summarize, one finally gets 

ttgwW * ^f^(k)T y 2 (k)T 2 (k)T^(k)V h {k) , (51) 

where the three transfer functions are given by Eqs. Q35p. 
(|4Tj). (|46| and (|50| and their respective pivot wavenum- 
bers by Eqs. K32J), ([21]), (gU) and @J. All of those four 
wavenumbers therefore correspond to frequencies of the 
gravitational wave / eq , / yon d, / xe nd and f rch . 

III. FORECASTS ON THERMAL HISTORY 
FROM FUTURE CMB AND GW EXPERIMENTS 

In this section, we discuss how and to what extent we 
can probe the thermal history of the Universe with future 
CMB and direct detection of GW experiments. 

A. An illustrative scenario 

For illustrative purposes, we consider in the follow- 
ing the specific scenario mentioned earlier. Inflation is 
driven by a massive field <p having the potential V{4>) = 
(l/2)m 2 2 with m being the inflaton mass. Furthermore 
there also exists another scalar field, denoted as a that 
comes to dominate the universe at later times. After in- 
flation, <j) oscillates around the minimum of its potential 
and the reheating (inflaton oscillation dominated) era is 
matter-like [27j till the Universe thermalizes and becomes 
radiation-dominated (the A-era). The reheating temper- 
ature T re h, associated with the energy density p re \i, refers 
to the time at which the reheating era ends and the ra- 
diation A-era starts. Then the field a begins to oscillate 
at some epoch and starts the a oscillation dominated 
era, referred to as the F-era in the previous section. We 
moreover assume that a decays after it dominates the 
Universe to start the usual radiation dominated period. 
Thus, the thermal history of the Universe proceeds as fol- 
lows: Inflation — > reheating era (oscillating (f> dominated 
era) — > RD (A-era) — > Oscillating a dominated era (Y- 
era) — > RD era. The final RD era continues until the 
radiation-matter equality, just before the recombination 
epoch. As discussed in the previous section, depending 



on the duration of each era, the predictions for CMB and 
GW spectra are different, from which we can probe the 
thermal history. 

Under these hypothesis, we can express all the quanti- 
ties of the previous section in terms of temperatures and 
frequencies. With R x — 1, assuming that g — g s = at 
the time of "yend," "xend," and "reh" and using Eq. (|22|) . 
Eqs. (JMI), (El and (gSJ), one obtains 

/ym ^ . 2Hz (^)(_^_), 

where T a is the temperature at which a decays into ra- 
diation. 

For CMB constraints, one can also simplifies the quan- 
tity i? ra d-Ry in terms of temperatures as 

ln(i? rad i? y ) = iln^-ilnF 

1 , Pond 1 , / 7T 2 ff*rchff*ycnd \ 

12 "A/* +12 11 {30 5 , xond )> 

(53) 

where u> rc h = has been used. Since T Ie h is assumed to 
be bigger than T a , we have the expected hierarchy / rc h > 
/xend > /yend- For the effective degrees of freedom, we 
assume those of the standard model of particle physics, 
i-e-, g* Ich = g* xcnd = 9*yend = 106.75 in the following 
analysis. Notice that this is certainly not verified but the 
influence of <?* remains very small in the final forecasts, 
the sensible quantities being logarithmic, see Eq. (|53l) . 

Some example plots for GW spectrum are shown in 
Fig. [3 The GW spectrum is obtained by numerically 
integrating Eq. (|21)|) . using the WKB approximation 
(h oc e ±lkv /a) for the oscillating phase. To see how the 
spectrum shape depends on the reheating temperatures 
and F, we show several cases for these parameters. For 
reference, the sensitivity frequency bands for BBO/FP- 
DECIGO and Ultimate DECIGO are also depicted (for 
the specifications of these experiments, see Table ITT|) . In 
some cases, the change of the spectrum can be traced 
with future experiments, in particular, for Ultimate DE- 
CIGO. In such a case, the parameters such as the reheat- 
ing temperatures and F can be well determined, which 
is going to be studied by using Fisher matrix analysis in 
the following. 

B. Fisher matrix analysis 

To forecast constraints from future experiments, we 
use a Fisher matrix analysis for both CMB (2|| and GW 
direct detection [2!|. Under the assumption of a Gaus- 
sian likelihood, the Fisher matrix is given by the second 
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FIG. 5: GW spectrum for several values of F ranging from 
F = to F = 100 at fixed T a (top). The bottom panel shows 
the spectra at fixed F = 10 and for various values of T CT . 
The reheating temperature T re h from the inflaton is fixed at 
T r eh = 10 12 GeV in these figures. 



derivative of the log-likelihood with respect to the pa- 
rameters pi at the likelihood maximum, 



d 2 InC 

dpidpj 



(54) 



and its inverse gives marginalized la error of the param- 
eter of interest, 



1. CMB experiment 



(55) 



In the case of CMB experiment, we assume that the 
Fisher matrix for both CMBpol and PLANCK is given 
by 



=2 XX' YY 



, dpi 



dp. 



(56) 



where X and X' are summed over the CMB temper- 
ature, E-mode polarization, and -B-mode polarization 
(X = T, E, B). The covariance matrix is given by 



EEEE — 



(Cov ^)tttt 

(Cov" 1 

(Cov~ 1 )bbbb 
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+ l)/sky 
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+ l)/sky 
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+ l)/sky 
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(21 


+ l)/sky 




Bi 2 )(Cl 
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(21 


+ l)/sky 
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(21 


+ l)/sky 
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(21 


+ l)/sky 



(Cf T 

(Pi 
(Pi 

[(CD 2 



-1 r>-2\2 
1 ^ ) 



tEE 



% 



{CJ E )\ 
CJ E (CJ T 



Cj E (Cf E 



w T B\ 



-2\2 



wT^BJ 2 



where " ' 1 



(T,P) 



= 4na 2 Tp) /N pix 



(57) 

is the variance of the 



noise temperature per pixel (in /J.K 2 ). For simplicity, we 
have assumed a Gaussian beam Bt ~ exp[— £(^ + l)cr 2 /2] 
with Ob = 0/ V8 In 2 being the beam width. In Table U 
we list the values of the observed fraction of the sky 
/sky, the temperature noise per pixel cry, the polariza- 
tion noise per pixel ap and the Gaussian beam width 
9 2 = 47r/7V pix for Planck [H and CMBpol [3l| exper- 
iments that are adopted to derive our forecasts. The 
different frequency channels are combined according to 
W( T P )i?| = wV (T.p)(B v t ) 2 ; where v refers to each 
channel component [32j. Finally, the maximum multi- 
pole value has been set to ^ max = 2000. 

Hereafter, we assume a flat ACDM Universe and set 
the fiducial cosmological parameters to the WMAP7 
mean values B . 



2. GW direct detection 

For the GW direct detection experiments, we will 
be considering three future experiments, FP-DECIGO, 
BBO and Ultimate-DECIGO FP-DECIGO is 

planned to be a Fabry-Perot Michelson interferometer, 
while BBO experiment will use time-delay interferometry 
(TDI) . Although they use different technology, in this pa- 
per, we do not distinguish these experiments, since they 
have similar sensitivit y. T he Fisher matrix for GW direct 
detection is given by [29( 



J~ij — 



3Hl 
10tt 2 



2 1 



bs 



(I,J) 



d/ 



\iu(f)\ 2 d Pi n gw (f)d P] n gw (f) 
rsj(f)s.j(f) 



(58) 
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TABLE I: Instrument parameter values for CMB experiments 

Experiment / s k y Center frequency [GHz] 9 [FWHM arcmin] ot [a*K] °"p [mK] 

Planck [30] 0.65 70 14 12.8 18.2 

100 10 6.8 10.9 

143 7.1 6.0 11.4 

217 5.0 13.1 26.7 

CMBpol [31] 0.65 100 4.2 0.87 1.18 

150 2.8 1.26 1.76 

220 1.9 1.84 2.60 



where f2 gw is given in Eq. (|24p and £ b s is the observation 
time. The subscripts / and J refer to independent sig- 
nals obtained at each detector, or observables generated 
by combining the detector signals. For a BBO-like ex- 
periment, the the summation runs over the TDI channel 
output index (/ = A,E,T) [29]. The overlap reduction 
Ju(f) for TDI variables is calculated following the pro- 
cedure of Rcf. [33]. For the cross term (I ^ J), we have 
"fu(f) = 0. The noise spectrum S7,j(/) is given by 



S E (f) = 8 sin' 
2 [3 + 2 cos / - 



; (//2){(2 
r- cos(2/) 



S T (/)=2(l + 2cos/ 



5' 



shot 



COS/)5 s hot 
ccclj' ) 

4sin 2 (//2)5 ; 



S : 



acccl 



(59) 

where / = 2irLf. The values of the arm length L, the 
shot noise S s hot and the radiation pressure noise Saccei 
have been reported in Table [H] for BBO/FP-DECIGO 
and Ultimate-DECIGO, respectively. One may be con- 
cerned about the noise contamination from white dwarf 
binaries at frequencies below ~ 0.1 Hz, and introduce a 
low-frequency cutoff to the integral. However, it may be 
possible to remove it by identifying all binaries and sub- 
tracting their contributions from data streams [3J-|36j. 
In this paper, we assume that this is the case and do not 
introduce such a low-frequency cutoff. 



C. Forecasts on early Universe history 

Now we investigate the constraint on the thermal his- 
tory of the Universe using the Fisher matrix method dis- 
cussed in the previous section. For this purpose, we focus 
on the reheating temperature T re h associated with the in- 
flaton, the temperature associated with the second field 
Tfj , and the amount of late-time entropy production F 
(or equivalently R y ). 



1. Unobservable spectral features 

In this section, we have performed a Fisher analysis to 
obtain the expected constraints from Planck/CMBpol for 



CMB and BBO/DECIGO for GW direct detection. As a 
fiducial model for the analysis, we consider two cases with 
9*){t T ™h ~ 10 16 GeV (corresponding to N = -AN* ~ 
57) and g^J^T^ ~ 10 9 GeV (corresponding to N ~ 52), 
assuming that there is no late-time entropy production, 
i.e., F=l. 

For the above fiducial models, there is no spectral sig- 
natures in the sensitivity range of the GWs direct de- 
tection experiments and the Fisher analysis remain in- 
sensitive to any information coming from the spectral 
shapes. As a result, one expects some degeneracies to 
occur between the model parameters. In Figs. [5] and 
we have represented the 2a allowed regions in the plane 
(T rc h,-F) and (m,F), respectively, for several combina- 
tions of those future data. 

There are two cases to consider according to the val- 
ues of T a . In Fig. [SI we have also represented the region 
where the signal-to-noise ratio is S/N > 5 for BBO/FP- 
DECIGO as the light blue shaded region. As studied 
in Fig. [5[ a late time entropy production would induce 
a suppression of the GW spectrum amplitude at fre- 
quencies higher than / ye nd, which corresponds to T a in 
Eq. (52). For BBO/FP-DECIGO, the suppression oc- 
curs in the sensitivity region when T a < 10 6 GeV. In 
that situation, a direct detection of the primordial GWs 
by BBO/FP-DECIGO would immediately yield a strong 
upper bound for F < 2.4 (see Fig. [5). At the same time, 
since the reheating also induces a suppression of the GWs 
spectrum amplitude, direct detection would readily ex- 
cludes T rc h < 10 6 GeV for the same reason. 

The opposite situation, i.e. the case with T a > 
10 6 GeV, the spectrum is suppressed at higher frequen- 
cies and beyond the sensitivity range of the experiment. 
Therefore, the amplitude of the observable GWs spec- 
trum remains mostly independent of the values of F. In 
this case, direct detection of GWs still allows a large pa- 
rameter space represented as the light yellow shading in 
Fig. [6l Within this region, the upper limit of F, equiv- 
alently the lower limit of T re h, is imposed by the con- 
sideration that reheating should end before the late-time 
entropy production begins. In other words, / XC nd < /rch 
gives T a F 2 / 3 < T^F- 1 / 3 using Eq. {52|), and given 
T a > 10 6 GeV, we get T leh /F > 10 6 GeV. Since the con- 
straints on F are completely degenerated with the one 
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TABLE II: Instrument parameters for GW experiments 
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FIG. 6: Future constraints on the (T re h,F) plane from 
CMB (Planck and CMBpol) and/or direct detection of 
GW (BBO/FP-DECIGO and Ultimate DECIGO). The light 
shaded region traces the 2a confidence interval of the two- 
dimensional marginalized posterior under the Fisher matrix 
analysis. For the darker blue shaded region traces the pa- 
rameter space in which the signal-to-noise ratio S/N > 5 (see 
text). 



bi? 



20 
18 
16 
14 
12 
10 
8 
6 
4 
2 




N=57 



/ BBO/FP-DECIGO 
Planck+BBO/FP-DECIGO 
/ ;? / CMBpol 




1.5 x 10 1 



2.0 x 10" 
m[GeV] 



2.5 x 10 1 



20 
18 
16 
14 
12 
10 
8 
6 
4 
2 



N=52 



BBO/FP-DECIGO 
Planck+BBO/FP-DEClGO 
CMBpol 

1 



1.5 x 10 



13 



13 



2.0 x 10 
m[GeV] 



2.5 x 10 1 



-2.1 
-0.1 

I. 9 

3.9 to 
5.9 3 

CD 

7.9 c5 

CD 

9.9 ^ 

II. 9 
13.9 
15.9 



-8.7 
-6.7 
-4.7 

-2.7 to 

-0.7 3 

a- 

1.3 
3.3 
5.3 
7.3 
9.3 



O 

tt) 
< 



FIG. 7: Future constraint in the plane (m, F) from CMB 
(Planck and CMBpol) and GW (BBO/FP-DECIGO), for the 
two fiducial models having T re h — 10 9 GeV (N = 52) and 
T rch ~ 10 16 GeV (N = 57). 



on T ro h, we have also reported the values of T re h along 
the right vertical axis of Fig. [71 

Regarding CMB constraints, there is also a degener- 
acy between T rc h and F, which is clear from Eq. ([551) . 
However, the infiaton mass m can be probed through the 
amplitude of primordial curvature perturbation while n s 



and r are related to F and T re h- On the other hand, 
the constraints from GW experiments basically come 
from their sensitivity to r and n T . Since in the single 
field inflation model, these two quantities are related by 
n T = — r/8 (at leading order in slow-roll parameters), 
m and F are constrained in a different way compared to 
CMB. In addition, CMB and GW direct detection exper- 
iments measure slow-roll parameters quantities at differ- 
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FIG. 8: Spectra of the three fiducial models (a), (b) and (c) 
together with the sensitivity region of the Ultimate-DECIGO 
GW experiment (see text). The corresponding forecasts are 
represented in Fig. [5] 



ent scales. As a result, the direction of this degeneracy 
differs thereby showing the complementarity of these ob- 
servables. 



2. Detection of spectral features 

For Ultimate-DECIGO, the suppression region of the 
GW signal occurs at T a < 10 3 GeV. However, there 
are large possibilities that the extreme sensitivity of 
Ultimate-DECIGO enables us to measure some of the 
spectral features. In that case, the detection of GWs 
would provide a precise determination of the parame- 
ters [38j . In the following, we investigate in more details 
the determination of the thermal history parameters by 
Ultimate-DECIGO in combination with CMBpol. 

In Fig. [9j we show the constraints on the parameters 
T rc h and F expected from Ultimate-DECIGO together 
with ones expected from CMBpol. Three different fidu- 
cial models have been considered: 

• (a) Both transition frequencies are inside the sen- 
sitivity frequency band of the experiment. We use 
the fiducial values T a = 10 4 - 9 GeV, T rch = 10 8 GeV 
and F = 10 19 . 

• (b) The frequency / xe nd is outside the sensitivity 
range. As such an example, we consider the case 
with T a = 10 3 GeV, T rch = 10 7 GeV and F = 10 2 . 

• (c) The frequency / rc h is outside the sensitivity 
range, which occurs for the fiducial values T a = 
10 4 - 5 GeV, T rch = 10 9 - 2 GeV and F = 10 17 . 

The spectra of these three models have been represented 
in Fig. [SJ Here, we take T re h, F, T a and m as free 



parameters and the predicted constraints are obtained 
by marginalizing over the remaining parameters. For 
CMBpol, the other cosmological parameters are also 
marginalized. 

In the case (a), all features of the reheating and late- 
time entropy production are within the GW sensitivity 
range. Since, in this case, Ultimate-DECIGO alone can 
well determine all the parameters, we do not show the 
constraint from CMBpol alone (in fact, the la contour 
from CMBpol ends up being outside the range of the fig- 
ure). However, it should be noted here that adding CMB 
data improves the constraint, as seen from the figure. 

For the case (b), the damping of the GW spectrum 
within the sensitivity zone is due to the inflationary re- 
heating and T ro h can be inferred. On the other hand, we 
find that GW measurements alone cannot determine the 
value of F, since the feature of late-time entropy produc- 
tion is now outside the region of detectability. However, if 
we combine CMB data, the constraint can be improved 
as seen from second panel of Fig. To illustrate this 
more clearly, we show in Fig. [TU] the marginalized con- 
straints for case (b) in the (m, F) plane. Constraints 
from direct detection have a strong parameter degener- 
acy between m and F, since they both cause suppres- 
sion of the amplitude of GWs at the maximum sensi- 
tivity region of Ultimate-DECIGO. However, CMB data 
such as CMBpol are strongly sensitive to the value of m 
and greatly helps to break the degeneracy. Therefore, al- 
though CMBpol data alone only provides constraints on 
F and T re h with large uncertainties (see middle panel of 
Fig. O, it can significantly tighten the constraints from 
Ultimate-DECIGO through improving the determination 
on m. 

In the case (c), the reheating frequency is outside the 
sensitivity region such that the value of T re h can no longer 
be well determined by GW measurements. Still T a and 
F can be well probed instead. As a result, the constraint 
from Ultimate DECIGO lies parallel to the T re h axis (see 
Fig.[S|). However, as repeatedly emphasized above, CMB 
probes the parameters differently from GW, thus in com- 
bination with CMBpol, T rc h can again be inferred and 
this ends up with breaking its degeneracy with F. Fur- 
thermore, in Fig. [TTJ we have also depicted the expected 
constraints in the (m,!^) plane, which also clearly illus- 
trates the complementarity of observations of CMB and 
GW. 



IV. SUMMARY 

We have discussed how one can probe the thermal his- 
tory of the Universe from the era just after the end of 
inflation until the BBN epoch. In any given inflation- 
ary models, the spectral index and the tensor-to-scalar 
ratio are related to the e-fold number at which a refer- 
ence scale exited the Hubble radius during inflation, thus 
in turn, they give a constraint on the total expansion of 
the scale factor since the above-mentioned epoch. By 
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FIG. 9: Predicted constraints in the (T Ic h,F) plane for the 
three different fiducial models (a), (b) and (c). The dotted 
and dashed lines show the marginalized 2a confidence con- 
tours expected from Ultimate-DECIGO and CMBpol, respec- 
tively. The solid line is the combined constraints. See also 
Fig. [TO] 
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assuming chaotic inflation, we have presented how the 
CMB is constraining the thermal history when there is 
an epoch of late-time entropy production, as for instance, 
a scenario where an oscillating scalar field dominates the 
Universe at some epoch, and then it decays. 

In the future, direct detection GW experiments are ex- 
pected to provide new cosmological probes, in addition 
to more precise CMB experiments such as Planck and 
CMBpol. Although CMB data can be used to determine 
the total amount of the cosmic expansion, they are only 
sensitive to the integrated thermal history and an epoch 
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of late-time entropy production remains fully degener- 
ated with the standard inflationary reheating era. This 
is not the case for GW experiments as they are precisely 
sensitive to the transition between these epochs. In par- 
ticular, for the scenario with late-time entropy produc- 
tion, the GW spectrum bends twice at some transition 
frequencies (see Figs. [5] and [5]). We have shown that if 
one, or more, frequencies do not lie in the GW experiment 
sensitivity region, CMB experiments can still greatly help 
to break the degeneracy. As shown in this paper, future 
experiments of GW and CMB are complementary in per- 
forming a tomography of the thermal history between the 
end of inflation and BBN. As such, they are expected to 
play a major role in our understanding of the whole his- 
tory of the Universe. 
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